Event-by-Event Analysis of Baryon- Strangeness Correlations: 
Pinning Down the Critical Temperature and Volume of QGP Formation 
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The recently proposed baryon-strangeness correlation (Cbs) is studied with a string-hadronic 
transport model (UrQMD) for various energies from _Ei a b = 4 AGeV to i/s = 200 AGeV. It is 
shown that rescattering among secondaries can not mimic the predicted correlation pattern expected 
for a Quark-Gluon-Plasma. However, we find a strong increase of the Cbs correlation function 
with decreasing collision energy both for pp and Au+Au/Pb+Pb reactions. For Au+Au reactions 
at the top RHIC energy (■s/s — 200 AGeV), the Cbs correlation is constant for all centralities 
and compatible with the pp result. With increasing width of the rapidity window, Cbs follows 
roughly the shape of the baryon rapidity distribution. We suggest to study the energy and centrality 
dependence of Cbs which allow to gain information on the onset of the deconfinement transition in 
temperature and volume. 
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Several observables 0] have been proposed throughout 
the last decades to study the characteristics of the highly 
excited matter created in heavy ions collisions, where a 
Quark-Gluon Plasma (QGP) is believed to be created. 
Among these observables, that give the opportunity to 
probe whether or not the system went through a phase 
of deconfined quarks and gluons, the ones related to fluc- 
tuations and correlations seem to be the most prospec- 
tive. Fluctuation probes might be more adequate for the 
exploration of heavy ion reactions, because the distri- 
butions of energy density or initial temperature, isospin 
and particle density have strong fluctuations from event 
to event 0, S 01 • O n the theoretical side event-by-event 
fluctuations where suggested to study 

• kinetic and chemical equilibration in nuclear colli- 

sions [EiasIliaHlEHH, 

• the onset of the deconfinement phase 0, 0, ITtI 

mm hi m 

• the location of the tri-critical end-point of the QCD 
phase transition m E3 or 

• the formation of exotic states, like DCCs [25j . 

On the experimental side, progress has been made by 
many experiments to extract momentum and particle 
number ratio fluctuations from heavy ion reaction: E- 
by-E fluctuations are actively studied in the SPS energy 
regime (starting from 20 AGeV on) by the NA49 group 

[MmmjliJMMpimill and thc CERES collab- 

oration J3| l3rflHl3a|~At RHIC energies thc PHENIX 
[H El El and STAR [H EH E| experiments are ad- 
dressing the field of single event physics. 

Recently a novel event-by-event observable has been 
introduced by Koch et al. the baryon-strangeness 

correlation coefficient Cbs- This correlation is proposed 
as a tool to specify the nature (ideal QGP or strongly cou- 
pled QGP or hadronic matter) of the highly compressed 
and heated matter created in heavy ions collisions. The 



idea is that depending on thc phase the system is in, 
the relation between baryon number and strangeness will 
be different: On the one hand, if one considers an ideal 
plasma of quarks and gluons, strangeness will be carried 
by freely moving strange and anti-strange quarks, car- 
rying baryon number in strict proportions. This leads 
to a strong correlation between the baryon number and 
strangeness. On the other hand, if the degrees of free- 
dom are of hadronic nature, this correlation is different, 
because it is possible to carry strangeness without baryon 
number, e.g. in mesons or QGP bound states. 

To quantify to which degree strangeness and baryon 
number are correlated, the following correlation coeffi- 
cient has been proposed 0] : 



Cbs 



, (BS)-(B)(S) 



(1) 



where B is the baryon charge and S is the strangeness. 
If a QGP is created, the expected value of Cbs will be 
unity as expected from lattice QCD, compatible with the 
ideal weakly coupled QGP. In the case of a hadron gas, 
where the correlation is non trivial, this quantity has been 
evaluated in to be Cbs = 0.66. 

In this paper, we study the correlation coefficient 
Cbs with the Ultra-relativistic Quantum Molecular Dy- 
namics model (UrQMD v2.2). The UrQMD is a non- 
equilibrium microscopic transport model that simulates 
the full space-time evolution of heavy ions collisions. It 
is valid from a few tens of McV to several TcV per nu- 
clcons in the laboratory frame. It describes the rescat- 
tering of incoming and produced particles, the excitation 
and fragmentation of color strings and the formation and 
decay of resonances. This model has been used before 
to study event-by-event fluctuations rather successfully 
EU HH US and yields a reasonable description of 
inclusive particle distributions. For a complete review of 
the model, the reader is referred to [4a,l47j . 
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Figure 1: Correlation coefficient for central Au+Au collisions 
at y/s = 200 shown as a function of the maximum rapidity ac- 
cepted. Circles are the calculation with HIJING [isll . Squares 
are the result of the UrQMD calculation and triangles of the 
JETSET for e+e~ at Vs = 200. The arrows are the values of 
a quark gluon plasma and of an hadron gas at a temperature 
T = 170 and chemical potential fi B = Hjj. Both HIJING 
and JETSET results are taken from 145^1 . 

Since the UrQMD is based on hadrons and strings 
it provides an estimate of the Cbs value in the case 
where no QGP is created, however taking into account 
the rescattering and the non-equilibrium nature of the 
heavy ion reactions. Cbs is evaluated from the event-by- 
event fluctuation analyses following [45t |48| : 

r A En B^SW - {± En En S™) 

^E„(S(")) 2 -(^E„^) 2 

(2) 

B^ n > and S^ n > stand for the baryon number and 
strangeness in a given event n. 

The correlation coefficient Cbs is depicted in Fig. 
as a function of the maximum rapidity accepted (\y\ < 
Umax)- The analyzed sample consists of central Au+Au 
events at y/s = 200 AGeV. For small acceptance windows 
around midrapidity, Cbs stays roughly constant. While 
for a large acceptance window, Cbs increases due to the 
inclusion of the fragmentation region with high baryon 
density. The different models deviate from each other for 
large acceptances due to differences in the handling of 
the fragmentation region, with small rapidity acceptance 
(relevant for the RHIC experiments), HIJING, JETSET 
and UrQMD yield consistent results. If the window ac- 
ceptance covers all produced particles, Cbs has to vanish 
because of baryon number conservation. 

In case a QGP is created, the signal given by the Cbs 
coefficient should survive the hadronic phase. With a 
strong enough longitudinal flow, strangeness and baryon 
number within a given rapidity range should be frozen in. 
The used rapidity window can not be too wide in order 
to avoid global baryon number and strangeness conserva- 
tion. Nevertheless, the acceptance window must be wide 



Figure 2: Correlation coefficient Cbs for central 
Au+Au/Pb+Pb (full symbols) and minimum bias p + p 
collisions (open symbols) as a function of y/s. The maximum 
rapidity accepted is y ma x = 0.5. 



enough to avoid smearing due to hadronization. 

Figure [2 depicts the energy excitation function of Cbs 
in both p+p and centrals Au+Au/Pb+Pb collisions. As 
discussed in |45j . Cbs increases with an increase of the 
baryon chemical potential (j,b when going to lower beam 
energies. With increasing collision energy, and therefore 
decreasing /is, Cbs g° cs down to Cbs ~ 0.4 at the 
highest RHIC energy available. Surprisingly, the general 
trend is the same for both p+p and Au+Au/Pb+Pb. 
Measuring the energy dependence of Cbs correlation 
around midrapidity might therefore allow to map out the 
onset of the QGP production. 

The dependence of Cbs on the number participants is 
studied in figure [3] The number of participants is de- 
termined via the scaled number of n~ 's in 4-7T geometry 
(N part = 0.2 < 7r~ >). This quantity is proportional to 
the overlap volume of the colliding nuclei and thus to the 
number of participants. The UrQMD model predicts a 
flat dependence of Cbs on centrality. Cbs ~ 0.4 from 
p + p to central Au + Au events. This is in strong con- 
trast with what is to be expected if the system enters 
a QGP phase at some centrality. In this case Cbs will 
increase (whether linearly or as a step function depends 
on the onset behaviour of the QGP phase) from periph- 
eral AA or pp towards central AA collisions. This might 
allow to extract in detail the volume dependence of the 
deconfinement transition at RHIC. 

To summarize, we have studied the dependence of the 
baryon-strangeness correlation coefficient as a function 
of the center of mass energy from = 4 AGeV to 
y/s = 200 AGeV for pp and central Au+Au/Pb+Pb re- 
actions. At yfs = 200 AGeV we have explored the cen- 
trality dependence of the Cbs correlation. Cbs is found 
to decrease from the lower energies towards the top RHIC 
energy available (here Cbs ~ 0.4). For minimum bias 
Au + Au events at y/s = 200, we predict a flat central- 
ity dependence of Cbs near midrapidity. At the highest 
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RHIC energy the Cbs value from the microscopic trans- 
port model is roughly half the one expected in the case 
of a QGP. We suggest to study the energy and centrality 
dependence of Cbs which allow to gain information on 
the onset of the deconfinement transition in temperature 
and volume. The CERES/NA49 and STAR experiments 
should be able to perform these analysis with their accu- 
mulated data. 



Figure 3: Correlation coefficient for Au+Au collisions at 
y/s — 200 as a function of the number of participants. The 
maximum rapidity accepted is Umax = 0.5. Full symbols are 
the results for Au+Au and the open symbol shows the p+p 
value. 
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